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Abstract 

By employing the QCD factorization approach, we calculate the new physics contributions to 
the branching radios of the two-body charmless B — > PV and B — > VV decays in the framework 
of the minimal supergravity (mSUGRA) model, we choose three typical sets of the mSUGRA 
input parameters in which the Wilson coefficient Ci~ f {mi ) ) can be either SM-like ( the case A and 
C) or has a hipped-sign (the case B). We found numerically that (a) the SUSY contributions 
are always very small for both case A and C; (b) for those tree-dominated decays, the SUSY 
contributions in case B are also very small; (c) for those QCD penguin-dominated decay modes, 
the SUSY contributions in case B can be significant, and can provide an enhancement about 
30% ~ 260% to the branching ratios of B — > K*(tt, (ft, p) and K<j> decays, but a reduction about 
30% ~ 80% to B — > K(p,u) decays; and (d) the large SUSY contributions in the case B may 
be masked by the large theoretical errors dominated by the uncertainty from our ignorance of 
calculating the annihilation contributions in the QCD factorization approach. 
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I. INTRODUCTION 



Along with the excellent performance of the B factory experiments 0, 0] and the forth- 
coming LHCb and other B meson related experiments 3], huge amount of B data with 
great precision are expected in the following five to ten years. The precision measurements 
of the B meson system can provide an insight into very high energy scales via the indirect 
loop effects of the new physics beyond the standard model (SM) |4|, |5j. Although cur- 
rently available data agree well with the SM predictions when considering the still large 
uncertainties both in theory and experiments, we generally believe that the B-factories 
can at least detect the first signals of new physics if it is there. 

As for the charmless hadronic two-body B decays considered here, people indeed have 
found some deviations from the SM expectations, such as the unexpected large branching 
ratios of B — > Krj' and B — > 7r°7r° decay modes, the " Kn" puzzle 5|, |(| and the so-called 
(f)K s anomaly at 3.8a level Although not convincing, these potential deviations may 
be considered as the first hints of new physics beyond the SM in B meson experiments. 

On the theory side, by employing the low energy effective Hamiltonian and various 
factorization hypothesis the charmless B meson decays have been widely studied in the 
framework of the SM |H 0, [l(l U, E 3- The possible new physics contributions 
to B meson decays induced by loop diagrams involving various new particles have been 
studied extensively, for example, in the Technicolor mo d els ll 511. the two-Higgs-doublet 
models [l^], 17] and the supersymmetric models [H, 1^, 2^, 2lL 221 . 

Among the various new physics models, the minimal super gravity (mSUGRA) model 
[23( | is a constrained minimal supersymmetric standard model [24| which has only five free 
parameters ( tan/3, mi , m , A , and sign(fi) ) at the high energy scale. In Refs. [3, [3, 

2ol . ijj . for example, the authors studied the rare decays B — > X s 7, B — > X s ll, B — > l + l~ 



and the B° — B° mixing in the mSUGRA model, and found some constraints on the 
parameter space of this model. 

In a recent paper [2(|, we calculated the SUSY contributions to the branching ratios 
of the twenty one B — > PP (P stands for pseudo-scalar light meson) decay modes in the 
mSUGRA model by employing the QCD factorization approach. We deduced analytically 
the contributions given by the new penguin diagrams induced by gluinos, charged-Higgs 
bosons, charginos and neutralinos, and obtained the analytical expressions of the SUSY 
corrections to the Wilson coefficients. By considering the current constraints, such as 
the branching ratio of the inclusive radiative B-decay, Cr 7 (mb) can be either SM-like or 
sign-flipped comparing with that in the SM. At the considered parameter point where the 
SUSY contributions can make the sign of CV 7 (rn&) reversed with respect to the SM one, 
we found that (a) the SUSY corrections to the Wilson coefficients C 7l (M w ) and C Sg (M w ) 
can be rather large; (b) the SUSY enhancements to those penguin-dominated decays can 
be as large as 30%-50%; (c)for B — > Kit decays, the inclusion of the SUSY enhancements 
can improve the consistency of the theoretical prediction with the data; and (d) the large 
SUSY contributions to B — > Krj' decays will help us to give a new physics interpretation 
for the so-called u Krj " puzzle. 

In this paper, we will extend the previous work |26[ to the cases for thirty nine B — > PV 
decays and nineteen B — > VV decay modes (here V stands for the light vector meson). 
This paper is organized as follows. In Sec. |H] and Sec. 11111 we give a brief review for 
the mSUGRA model and the QCD factorization approach for two-body B meson decays. 
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In Sec. HXl and Sec. El we present the related formulae and calculate the CP-averaged 
branching ratios for B — > PV and B — > VV decays, respectively. The summary and some 
discussions are included in the final section. 



II. THE MSUGRA MODEL AND THE NEW PHYSICS CONTRIBUTIONS 

In this section, we first recapitulate the basic theoretical framework of the mSUGRA 
model to set up the notation and then calculate the SUSY corrections to the Wilson 
coefficients in the mSUGRA model. 



Outline of the mSUGRA model 



The most general superpotential compatible with gauge invariance, renormalizability 
and R-parity conserving in MSSM can be written as 24 \: 



W 



[I) 



Where Q, U, D, L, E, Hi and H 2 are chiral superfields. fo, fu and Je are Yukawa coupling 
constants for down-type, up-type quarks, and leptons, respectively. The suffixes a, (3 = 
1, 2 are SU(2) indices and i,j=l,2,3 are generation indices. e a p is the antisymmetric tensor 

with 8x2—1. 

In addition to the SUSY invariant terms, a set of terms which explicitly but softly 
break SUSY should be added to the supersymmetric Lagrangian. A general form of the 
soft SUSY-breaking terms is given as: 



- C soft = (mg).. q^qLj + (mDij u* Ri u Rj + (m 2 D ).. d* Ri d Rj + (m 2 L ).. l + J Lj 
+ (ml)^ e* Ri e Rj + Ajh+h + A 2 2 h+h 2 



Auij&lu*^ + Ao^ql^ + A Eij h<fl p Li e* Rj + B^h p 2 

1 ,z„z. 1 

2 



+^m^BB + ^m w WW + ^m & GG H.C '. 



(2) 



where qu, u* Ri1 d* Ri , In, e* Ri , and hi and h% are scalar components of chiral superfields 
Qi, Ui, Di, Li, Ei, Hi, and H 2 respectively, and B, W, and G are U(l)y, SU(2)l, and 
SU(3)c gauge fermions. And the terms appeared in Eq.(J2J) are the mass terms for the 
scalar fermions, mass and bilinear terms for the Higgs bosons, trilinear coupling terms 
between sfermions and Higgs bosons, and mass terms for the gluinos, Winos and binos, 
respectively. 

In the mSUGRA model, a set of assumptions are added to the MSSM. One underlying 
assumption is that SUSY-breaking occurs in a hidden sector which communicates with the 
visible sector only through gravitational interactions. The free parameters in the MSSM 
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are assumed to obey a set of boundary conditions at the Plank or GUT scale: 



Oil — a 2 = a 3 = a x , 
Aj = A 2 2 = ml 

AjJij — fllijAo, Aoij = foijAo, Afiij = f'EijAo, 

nig = m^y = itIq = mi (3) 

where a, = gf/^An), while g^ (i=l,2,3) denotes the coupling constant of the U(l)y, 
SU(2)l, SU(3)c gauge group, respectively. The unification of them is verified according 
to the experimental results from LEPljl^] and can be fixed at the Grand Unification Scale 
Mqut ~ 2 x 10 16 Gev. 

Besides the three parameters mi, m and A , the bilinear coupling B and the supersym- 
metric Higgs(ino) mass parameter \x in the super symmetric sector should be determined. 
By requiring the radiative electroweak symmetry-breaking (EWSB) takes place at the low 
energy scale, both of them are obtained except for the sign of /i. At this stage, only four 
continuous free parameters and an unknown sign are left in the mSUGRA model |28| |: 

tan P, mi , mo, Ao, sign(fi) (4) 

In present analysis, we assume that they are all real parameters. Therefore there are no 
new CP-violating complex phase introduced other than that in the Cabbibo-Kabayashi- 
Maskawa(CKM) quark mixing matrix |2jj]. Once the five free parameters of the mSUGRA 
model are determined, all other parameters at the electroweak scale are then obtained 
through the GUT universality condition and the renormalization group equation (RGE) 
evolution. Like in Ref. j2^|, we here also use the Fortran code - Suspect version 2.1 |3p| - 
to calculate the SUSY and Higgs particle spectrum 1 . 



B. New physics effects in the mSUGRA model 

As is well-known, the low energy effective Hamiltonian for the quark level three-body 
decay b — > qq q (q G {d, s}, q G {u, c, d, s, b}) at the scale jtt~mj reads 

r f 2 



i=i 

10 



-v tb v; q CjtoOM - v tb v* {c 7 mo 7 m + C 8g (n)0 8g (n)] \ (5) 

j=3 



where V p bV* q with q — d, s is the products of elements of the CKM matrix [29( . The 
definitions and the explicit expressions of the operators Oj (i = 1 ~ 10, 77, 8g) and the 
corresponding Wilson coefficients Cj can be found in Ref. |3l| . In the SM, the Wilson 



1 For more details, one can see discussions in Ref. |3| and references therein. 
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coefficients appeared in eq.© are currently known at next-to-leading order (NLO) and 
can be found easily in Ref. [3l| . 

In the mSUGRA model, the new physics effects on the rare B meson decays will 
manifest themselves through two channels. 

• The new physics contributions to the Wilson coefficients of the same operators 
involved in the SM calculation; 

• The other is the new physics contributions to the Wilson coefficients of the new 
operators such as the operators with opposite chiralities with 0, appeared in eq.flSJ). 

In the SM, the latter is absent because they are suppressed by the ratio m s /mb. In the 
mSUGRA model, they can also be negligible, as shown in Ref. [32]. Therefore we here 
use the same operator base as in the SM. 

It is well known that there is no SUSY contributions to the Wilson coefficients at the 
tree level. At the one-loop level, there are four kinds of SUSY contributions to the quark 
level decay process b — > qq q , depending on specific particles propagated in the loops: 

(i) the charged Higgs boson H ± and up-type quarks u, c, t; 

(ii) the charginos X12 an d the scalar up-type quarks u, c, t; 

(iii) the neutralinos Xi 2 3 4 an d t ne scalar down-type quarks d, s, b; 

(iv) the guinos g and the scalar down-type quarks d, s, b. 



In Ref. [26], we have given a detailed derivation of the lengthy expressions of SUSY 
corrections induced by those penguin diagrams involving SUSY particles to the Wilson 
coefficients. 

In the mSUGRA model, there are five free parameters: tan /3, mi , m , A and sign(fi). 
Before we calculate the new physics contributions to the considered B — > PV, VV de- 
cays, we have to check the experimental bounds on these parameters, which have been 
studied extensively by many authors j33[ 13 I^H |36|. The supersymmetry parameter 
analysis(SPA) project with the main target of high-precision determination of the su- 
persymmetry Lagrangian parameters at the electroweak scale, for example, is under way 
now 35] . By confronting the precision electroweak data provided by the LEP experiments 
with the theoretical predictions within the mSUGRA model, the strong constraints on 
the parameter space have been obtained |33|, 34], 35l. 136^. 



From the well measured B — > X s 7 decays, the magnitude-but not the sign- of the Wil- 
son coefficient CV 7 (m&) is strongly constrained. This is an important constraint. However, 
various new physics can change the sign of CV 7 (mb) without changing the branching ratio 
of B — > X s 7 decay obviously. 

From the measurements of the b — > sl + l~ decays, one can determine the relative sign 
of the Wilson coefficients as well as their absolute values. The latest Belle and BaBar 
measurements of the inclusive B — > X s l + l~ branching ratios indicated that the sign of 
CVyfmft) is unlikely to be different from that in the SM [37]. On the other hand, the studies 



37t l38| also show that a positive could be made compatible with experiments only 

by large 0(1) new physics corrections to Cg^/ , while the SM values of Cg^/ are around 
+4.2 and —4.4 respectively. 
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TABLE I: Three typical sets of SUSY parameters to be used in numerical calculation. The 
third case is the "typical" mSUGRA point "SPS lb" as defined in Ref. [3^ . The last column 
show the values of the ratio R 7 = C 7 y(mb) / CyJ^ 1 (m,b) . All masses are in units of GeV. 



CASE 


uiq mi 


A 


tan/3 Sign[fi] 


R 7 


A 


300 300 





2 


1.10 


B 


369 150 


-400 


40 + 


-0.93 


C 


200 400 





30 + 


0.82 



For the semi-leptonic Wilson coefficients Cg and C\g in the minimal supergravity 
model, the authors of Ref. have made a detailed analysis and found that Cg{mb) 
and Cio{rrib) differ from their SM values by at most 5% in the parameter space for 
—2 < C 7l {mb) /Cj^' 1 \m h ) < 2 and for tan/3 = 3,30. This result supports a SM-like 
C 7l {m h ) in the mSUGRA model. 

This August, Belle Collaboration reported their measurement of the ratios of Wilson 
coefficients Ag/A 7 and A 10 /A 7 2 in B — > K*l + l~ decay They exclude a positive 

at more than 95% CL., but can not determine the sign of A 7 (i.e. C 7l (jrib))- 

Based on a lot of previous studies about the constraints on the parameter space of 
the mSUGRA model jiol I33L l34l I35I . l3ol l38l | , we can choose several typical sets of input 
parameters to show the pattern of the new physics corrections to the branching ratios 
of the studied decays in the mSUGRA model. In Ref. |26[ we selected two sets of input 
parameters (case-A and case-B) and found that (a) the SUSY corrections to the Wilson 
coefficients C^(k = 3 ~ 6) are always very small and can be neglected safely; (b) with 
the inclusion of SUSY corrections, the Wilson coefficient C 7l (rrib) remained SM-like for 
case-A, but changed its sign for the case-B. Numerical results for the branching ratios of 
B — > PP decays also show that the new physics corrections are very small for the case-A, 
but can be significant for the case-B. 

In this paper, besides the choice of case-A and B as defined in Ref. |26| , we also consider 
the third case, the case-C, which is the "typical" mSUGRA point "SPS lb" as defined in 
Ref. [36l| . In the case-C, the SUSY corrections to C^ik = 3 ~ 6) are also negligibly small, 
while CV 7 (m.&) is SM-like after the inclusion of SUSY contributions. The three sets of 
mSUGRA input parameters to be used in numerical calculations here are listed in Table 
HI the numerical values of the ratio R 7 = C 7l {m h ) / C 7 ^ (nib) for each cases are also given 
in the table. It is easy to see that the Wilson coefficient C 7l in the mSUGRA model 
is SM-like (negative) for both case-A and C, but nonstandard (positive) for case-B. The 



2 At next-next-to-leading order (NNLO), the effective Wilson coefficients ,Cg and C^q have 
many small high order correction terms 39], one usually use the leading coefficients A 7 , Ag and Aiq in 
the evaluations 01 . The numerical differences between C^g 10 and A 7t g t io are indeed very small. It 
is worth noting that, instead of the (electroweak penguin) coefficients Cg and C10 as defined in Eq.JSJ), 
the coefficients Cg^ and C^g here correspond to the low-energy interaction terms (sz^&l)^^) and 
(sl7 a1 6l)(/~7 ai 7 5 0, respectively. 
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ratio i?7 > 1 for case-A, but < 1 for the case-C, due to the strong bound imposed by the 
precision data of B — > X s ^. 

From above discussions, one expects that the case-A and C will be very similar phe- 
nomenologically. The numerical calculations to be given in the following sections indeed 
show that the theoretical predictions for the branching ratios in case-A and C are almost 
identical with each other: the difference is less than 2%. We therefore will present the nu- 
merical results of the branching ratios explicitly for case-B and case-C only, and compare 
them with the corresponding SM predictions. 



III. B M X M 2 DECAYS IN THE QCD FACTORIZATION APPROACH 

In QCD factorization approach, when the final state hadrons of B meson two-body 
decays are all light mesons (mesons composed of light u, d, or s quarks, and with a mass 
of order Aqcd), the matrix element of each operator in the effective Hamiltonian 7i e // 
can be written as he 



{M 1 M 2 \O i \B) = Vif^ Ml / dxT^{x)^ M M + {Mi M 2 ) 

+ f d£ / dx I dyT{ I ^x 1 y)^ B {i)^M 1 {x)^ M M (6) 
Jo Jo Jo 

where p^Mi - g ^ e f orm factor describing B — > M\ decays, and T^ 1 denote the per- 
turbative short-distance interactions and can be calculated by the perturbation approach, 
and §x{z){X = B,Mip) are the universal and nonperturbative light-cone distribution 
amplitudes (LCDA) for the heavy B meson and the light meson respectively. Weak 
annihilation effects which are suppressed by Aqcd / 'wj are not included in Eq.©. 

Considering the low energy effective Hamiltonian Eq.fjSJ) and the QCDF formula Eq.©, 
the decay amplitude can be written as 

A f (B - M 1 M 2 ) = W^WI^I 5 )^ ( ? ) 



p=u,c 



Here {MiM<2,\Oi\B) p is the factorized matrix element. The explicit expressions for the 
decay amplitudes of B — ► M\M 2 decays can be found for example in the Appendixes of 
Refs. HQ. Following Beneke et al. Q, the coefficients a i {M 1 M 2 ){i = 1 to 10) in Eq.© 
with Mi absorbing the spectator quark is 



a\ [M X M 2 ) = (d + ^)N i {M 2 ) 



47T 



2 



Vi(M 2 ) + — H i (M 1 M 2 ) 



2), 



where the upper (lower) signs apply when i is odd (even). The functions V^(M 2 ) account 
for one loop vertex corrections, Hi(M\M 2 ) for hard spectator interactions and Pf{M 2 ) 
for p eng uin contributions. The explicit expressions for these functions can be found in 
Rcf. 131. 
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As mentioned previously, the SUSY contributions to the Wilson coefficients of the 
four-quark penguin operators are very small and have been neglected. The new magnetic 
penguin contributions in the mSUGRA model can manifest themselves as radiative cor- 
rections to the coefficients etf and be contained in the functions P[(M 2 ), which is present 
only for i = 4,6, 8, 10 To make this visible, here we just show the functions Pf (M 2 ) 
and Pi! 



10 (M 2 



At order a s , these two functions are 



P£(M 2 



Cpa s 
+C 3 



4 ! m 6 2 

-In — + - - Gm 2 {s p ) 

6 /I 6 



6/2 6 



+(C 4 + c 6 



An 



n. 



2)G m M-Gm 2 {sc)-G M2 {1) 



dx 
1 — x 



(9) 



a: 



9irN r 



(Ct + N C C 2 ) 



'4 m b 2 

3ln- + --G M2 ( Sp ) 



dx 



x 



(10) 



where s u = m^/ml « and s c = m 2 c /m1 are mass ratios involved in the evaluation 
of penguin diagrams. &m 2 { x ) is the leading twist LCDA which can be expanded in 
Gegenbauer polynomials. C^l = CV 7 — |C 5 — C 6 and = C 8g + C 5 are the so-called 
"effective " Wilson coefficients, where the SUSY contribution is involved. The explicit 
expressions of the functions Gm 2 (0), Gm 2 (X) an d Gm 2 (s p ) appeared in Eqs. (|9ll0|) can be 
found easily in Refs. [ll|, [l3| . 

When calculating the decay amplitudes, the coefficients ai{i = 3 ~ 10) always appear 
in pairs. So in terms of the coefficients af , one can define af as follows [l3j : 



a x {M x M 2 ) 
a 2 {M x M 2 ) 

a p 3 (M x M 2 ) 



oP 4 (M x M 2 ) 



a 



3,EW 



a 



4,EW 



{M X M 2 



(M 1 M 2 



a x (M x M 2 ), 
a 2 (M x M 2 ), 

f af (MiM 2 ) - og(MiM 2 ) 
I a|(M!M 2 ) + af (MxM 2 ) 

<(M X M 2 ) + rf ^ ^(MiMs; 
<(MxM 2 )-rf ^(AfiM 2 ; 



if .\/ : .U> 
if M X M 2 



PP, VP . 
PV, 



if M X M 2 
if MiM 2 



PP, PV . 
VP, 



a%(M x M 2 ) - a p 7 (M x M 2 ) 
ag(MiM 2 ) + a p 7 (M x M 2 ) 

a? (M 1 M 2 )-rf 2 af(M 1 M 2 ) 



if MiM 2 
if MxM 2 

if MxM 2 
if MiM 2 



PP, VP 
PV, 



PP, PV 
VP. 



11) 
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For pseudoscalar meson P and vector meson V, the ratios and r^(//) are defined as 

r " = m 6 ( / i)(m, 1 +m„)( / i)' ( ' 

r v w = 2n»fM, (13 ) 

where m qi and m q2 are the current masses of the component quarks of P meson, and 
fy(fi) is the scale-dependent transverse decay constant of vector meson V. Although all 
the terms proportional to r^ 2 are formally suppressed by one power of Aqcd /nib in the 
heavy-quark limit, these terms are chirality enhanced and not always small. They are very 
important in those penguin- dominant B meson decays, such as the interesting channels 
B — > Kr] ,etc. 

In QCD factorization approach, the nonfactorizable power-suppressed contributions are 
neglected. However, the hard-scattering spectator interactions and annihilation diagrams 
can not be neglected because of the chiral enhancement. Since they give rise to infrared 
endpoint singularities when computed perturbatively, they can only be estimated in a 
model dependent way and with a large uncertainty. In Refs. 11, 3 these contributions 
are parameterized by two complex quantities, %h and xa, 

XH,A=^+PH,A^ A )^ T ^ (14) 

where A h = 0.5GeV, 4>h,a are free phases in the range [—180°, 180°], and Ph,a are real 
parameters varying within [0, 1]. In this paper, we use the same method as in Refs. mm 
to estimate these two kinds of contributions. 

As given in Ref. 0, the annihilation amplitude can be written as 



A ann {B ^ ^ <fF £ £ V^J B j M J u^M^) (15) 



p=u,c 



where fs and fu are the decay constants of B meson and final-state hadrons respectively. 
The coefficients bi{MiM 2 ) describe the annihilation contributions. For explicit expressions 
of coefficients b iy one can see Refs. (lli fl3| . 

For B — > VV decays we have two additional remarks: (a) since < V|gig 2 |0 >= 0, 
B — > VV decays do not receive factorizable contribution from and as penguin terms 
except for spacelike penguin diagrams; and (b) unlike the PP and PV decay modes, the 
annihilation amplitude in the VV mode does not have a chiral enhancement of order 
Mg/(m q mb). Therefore, it is truly power suppressed in heavy quark limit and will be 
neglected in our calculation. The explicit factorizable coefficients Oj for B — > VV can be 
found in Ref. El • 



IV. BRANCHING RATIOS OF B — > PV DECAYS 

In the following two sections we will calculate the CP-averaged branching ratios for 
thirty nine B — > PV and nineteen B — > VV decay modes, respectively. We usually use 
the central values of the input parameters as collected in Appendix A, and consider the 
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effects of the uncertainties of these input parameters as specified in the text for individual 
decay channels. 

Theoretically, the branching ratios of charmless decays B^PV in the B meson rest 
frame can be written as 

Br(B^PV) = ^\^\A(B^PV)\ 2 , (16) 
87T m B 

where t# is the B meson lifetimes, and \P C \ is the absolute values of final-state hadrons' 
momentum in the B rest frame and written as 

ip I _ Vl m B ~ ( m P + m v ) 2 }[m 2 B - (m P - m v ) 2 } 

1 c| 2m B ' 1 ' 

For the CP-conjugated decay modes, the branching ratios can be obtained by replacing 
CKM factors with their complex conjugate in the expressions of decay amplitudes. Using 
the decay amplitudes as given in Refs. @, El and the coefficients a, in Eq. ([5|) or a.i in 
Eq.()lljl. it is straightforward to calculate the CP-averaged branching ratios of those thirty 
nine B — > PV decay modes in the SM and mSUGRA model. 

From Eqs.(|5HTT|l. we can find that the potential SUSY contributions are mainly embod- 
ied in a^(MiM 2 ) and a^ ew (MiM 2 ). Therefore, one naturally expect a large new physics 
corrections to those penguin dominated B meson decays. 



Numerical results 



In Table HTland IIIIl we show the theoretical predictions for the CP- averaged branching 
ratios for B — > PV decays in both the SM and mSUGRA model, assuming p, = m&/2, rrib 
and 2rrib, respectively. And in the SM we give both Br^ +a and Br? , the SM predictions 
for the branching ratios with or without the inclusion of annihilation contributions, re- 
spectively. In the mSUGRA model, we consider both case B and case C and only give Br* 
so that we can compare the relative size of the new physics contribution with the annihila- 
tion contribution for each considered decay channel. From the numerical results, one can 
see that the SUSY corrections to the b — > s transition processes are generally larger than 
those to the b — > d processes because of the CKM factor suppression (|Vt&V^| ~ 10~ 2 ) in 
b —>■ d penguin transition. 

We classify the thirty nine B — ► PV channels into b —>■ d and b — > s processes. For 
b — » d processes, we have the following remarks 

• The B° — > 7r ± p =F and B~ — ► -K p~,-K~p ,-K~LU,ri^p~ decays. 

These channels are tree-dominated decay modes and depend on the large coefficient 
a\. The SUSY corrections to these decays are very small and can be neglected 
safely in all the parameter space. 

• The decays B° — > n°p°, tt°cu, r]^p°, rf^u). 

These channels have small branching ratios. The reason for this is twofold. On the 
one hand, the tree contributions of these channels are involved in the coefficient a 2 ~ 
0.2, which is far smaller than the large coefficient a.\ ~ 1. The penguin contributions, 
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TABLE II: Numerical predictions in the SM and mSUGRA model for CP-averaged branching 
ratios ( in units of 10~ 6 ) for b — » d transition processes of B — ► PV decays, where Br^ +a and 
Br$ denotes the branching ratios with and without the annihilation contributions respectively. 
For Case B and C, the branching ratios without annihilation contributions are given. 



D 

£3 - 


— > rv 


p = m b /2 




1* = 


m b 




p = 2mt 




-> dj 


SM 


mSUGRA 


SM 


mSUGRA 


SM 


mSUGRA 






Br 1 




(B) 


(C) 


Br 1 




(B) 


(C) 


B r J 


7-» f-X-fi 
J3y>J^T a 


(B) (C) 


B~ — 


_ n 

> 7T p U 


11.2 


11.2 


11.2 


11.2 


11.5 


11.4 


11.5 


11.5 


11.8 
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on the other hand, are strongly suppressed by the CKM factor |Vt&V^| ~ 1CT 2 . 
From Table |Hj we can see that the SUSY corrections are always smaller than the 
annihilation contributions, and thus can be easily masked by it. 

• The decays B° — ► 7r°0, rf^cj) and B~ — > n~(f). 

This kind of channels have both penguin and weak annihilation contributions. But 
the penguin contributions come from the small coefficients and a^ ew . Therefore 
weak annihilation contributions are dominant for B° — > r/^cf). The branching ratios 
of these decay are at the O(10~ 9 ) level, and the SUSY corrections in all parameter 
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TABLE III: Numerical predictions in the SM and mSUGRA model for CP-averaged branching 
ratios (in units of 10~ 6 ) for b — > s transition process of B — ► PV decays, where Brf +a and Br? 
denotes the branching ratios with and without the annihilation contributions respectively. For 
Case B and C, the branching ratios without annihilation contributions are given. 
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space can hardly affect them. 

• The decays B° -> K°K*° and B~ -> K~K*°. 

These channels are penguin dominant decays. In their amplitudes, the dominant 
term is proportional to af (P, V) and af ew (P : V) . The SUSY corrections therefore 
can increase their branching ratios significantly in some parameter space, such as in 
case B where the SUSY contributions can provide a 130% enhancement and are far 
larger than the annihilation contributions. But in case C where C 7l (m b ) is SM-like, 
the SUSY contributions are negligibly small. 

• The decays B° -> K°K*° and E>- -> K°K*-. 

Different from the B° — > K°K*° and B~ — > K~K*° decays, the dominant terms 
here are proportional to ccf (V, P) and a% ew (V,P). The SUSY corrections interfere 
destructively with their SM counterparts and will decrease the branching ratios by 
about 50% in case B. For case C, the SUSY contributions are still very small. 

• The decays B° -> K + K*~,K-R* + . 
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These two channels have weak annihilation contributions only. Therefore in the 
mSUGRA model their branching ratios can hardly be affected. 

For b —>■ s transition processes, the tree contribution is suppressed by CKM factors, 
the penguin contributions therefore play the major role. These decays can be classified 
into three groups. And our remarks are 

• B —> irK*, K(p decays. 

These decays are penguin dominant decays and both the annihilation contributions 
and the SUSY contributions can give enhancements to their branching ratios. But 
for case C the SUSY contributions are smaller than the annihilation contributions 
and can be masked easily. Only in case B where the enhancements can reach as 
large as 100% ~ 260% for B — > nK* decay and about 100% for B — > Kef) decay. 

• B — > Kp, Kuj decays. 

In case B the branching ratios of these decays will be decreased by 30% ~ 60% after 
the inclusion of SUSY corrections and theoretically we can isolate the SUSY contri- 
butions from the annihilation contributions for which tend to increase the branching 
ratios greatly. In case C, however, the SUSY contributions are small, while the an- 
nihilation contributions are relatively important for explaining the current large 
measured data. 

• B — > K*rj"> decays. 

For B — > K*rj decay, their amplitudes strongly depend on af (P, V) and af ew (P, V). 
Although the inclusion of SUSY corrections will increase their branching ratios, 
the SUSY contributions will also be masked by large annihilation contributions. 
For B — > K*7] decay, however, their amplitudes mainly depend on af (V, P) and 
a± ew (V,P), the SUSY corrections will trend to decrease their branching ratios in 
both cases and can be separated from the annihilation contributions which can 
contribute a 40% ~ 60% enhancement to the branching ratios. 

From above discussions, one can see that the SUSY contributions in case-C (preferred 
by the measured branching ratio of B — > X s l + l~ decay are always negligibly small. 



B. The data and phenomenological analysis 

Among the thirty nine B — > PV decay modes considered here, eighteen of them have 
been measured experimentally. The latest individual measurements as reported by dif- 
ferent groups and the new world average for the branching ratios can be found in the 
HFAG (Heavy Flavor Averaging Group) homepage 4^1 • In this subsection we will make 



a phenomenological discussion for those eighteen measured decay channels. 



1. B — ► up and B — * ttuj 



Among the B — > up and tcuj decays, four of them have been well measured, and the 
experimental upper limits are available for the remaining B — > 7r°p° and 7r u decay modes. 
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The data and the theoretical predictions for the four measured decay modes (in units of 
10" 6 ) in the SM and mSUGRA model (both Case B and Case C) are 
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Data, 
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(19) 



(20) 



(21) 



where the major errors are induced by the uncertainties of the following input parameters: 



B- 



0.28 ± 0.05, A 



0.37 ± 0.06. 



0.33 ± 0.05 and 7 = 57.8° ± 20 c 



Throughout this paper we take the central values and the ranges of these input parameters 
specified here as the default values, unless explicitly stated otherwise. 

We also set pa = as the default input value |llj in numerical calculations, and scan 
over pa £ [0, 1] and 4>a £ [—180°, 180°] to estimate the theoretical error induced by the 
uncertainty of annihilation contribution. 

For the convenience of analysis, in Fig. ^ we show the 7 dependence of the theoret- 
ical predictions 3 for the branching ratios of the three B — > np and B~ — ► h~oj decays 
respectively and the experimental data are also marked. 

From the numerical results as given in Eqs. (|18H21j) and Fig. Q one can see that 

• For these tree-dominated decay modes, the SUSY corrections in the mSUGRA 
model are always very small for 7 ~ 60°. 

• As to the theoretical errors, generally the uncertainties of the form factors and the 
annihilation contribution are the dominant sources. For B° — ► 7i ± p T , the annihila- 
tion contribution which has been parameterized by xa also give large uncertainty. 
For B~ — ► 7r~p°, however, the uncertainty of the annihilation contribution is small 
and has been ignored here. 



3 The central values of all input parameters except for the CKM angle 7 are used in this and other 
similar figures. The theoretical uncertainties are not shown in all such kinds of figures. 
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FIG. 1: The 7 dependence of the branching ratios of B — ► up and B~ — ► 7r _ w decays in the 
SM and mSUGRA model. The solid, dots lines and the scatter of prong show the central values 
of the SM prediction and the mSUGRA ones in both case B and case C, respectively. The 
horizontal back-slashed gray bands show the data. 

• The theoretical predictions for the branching ratios in both the SM and mSUGRA 
model are all consistent with the data within one standard deviation since the the 
theoretical errors are still large. 



2. ttK*, B^K4> 

Like B — > irK, the four B — > nK* decays are penguin dominated decay modes and 
therefore sensitive to the new physics contributions. The CP-averaged branching ratios 
of the four B — > nK* decays have been measured experimentally and the data are much 
larger than the SM predictions. In the mSUGRA model, the new penguin diagrams may 
contribute effectively to these decays. The data and the numerical results (in units of 
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FIG. 2: The same as FigJT]but for B — ► txK* decays. 
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where the renormalization scale \i varies from m^/2 to 2m;,, and the second error in Eq.( 
and (|23j) is induced by the uncertainty of the form factor A B ^ K , A B ^ K * = 0.45 ± 0.07. 

Fig. El shows the 7 dependence of the branching ratios for the measured B — > ttK* 
decays in both the SM and the mSUGRA model. From this figure and Eqs. fl22M23)l . one 
can see that 

• The central values of the SM predictions for the branching ratios are only about 
20 to 30 percent of the measured values. The SUSY contributions are also very 
small in the parameter space where CV 7 (mft) is SM-like. But for the case B, the 
SUSY contributions are large and can provide a factor of two enhancements to these 
penguin-dominated decays. After the inclusion of the large SUSY contributions, the 
mSUGRA predictions for B — > tt~K*°, tt°K*~ and tt°K*° become consistent with 
the data within one standard deviation. 

• For all the four B — > ttK* decays, the dominant error is induced by the uncertainty of 
the annihilation contribution, it is large and may mask the large SUSY contributions 
in Case B. 

• After the inclusion of the large theoretical errors, the theoretical predictions for the 
branching ratios of B — > tt~K*°, tt°K*~ and ir°K*° in both models can be consistent 
with the corresponding data within one standard deviation. But for B° — > tt + K*~ 
decay, only when we include the large SUSY contribution in Case B as well can the 
theoretical value match with the data within one standard deviation. 

Like B — > ttK*, the SUSY contributions to B — > Kef) in case B are large and can provide 
about 100% enhancement. From Fig. El one can see that the mSUGRA predictions in 
case B can become consistent with the data within one standard deviation. The data, the 
central values and the major errors of the branching ratios (in units of 10~ 6 ) in the SM 
and mSUGRA model are 
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FIG. 3: The same as Fig^but for B — ► Kcj) decays. 



From these numerical results, one can see that (a) the central values of the SM predic- 
tions are only about one third of the measured values. But the mSUGRA predictions in 
case B can become consistent with the data within one standard deviation, (b) Similar 
to the B — > 7rK* decays, the dominant error here also comes from the uncertainty of the 
annihilation contribution and the form factor F^ K . The error induced by xa are very 
large and even mask the large SUSY contributions in Case B. 



3. B —> Kp, B — > Klo 

Among the four B — ► Kp decays, three of them have been well measured experi- 
mentally. The data, the theoretical predictions and the major errors (in units of 10 -6 ) 
are 



Br(B° -> K~p* 



Br(B~ -> R-p ^ 



Br(B° -> K°p°" 



+1.6 
1.5' 



9.9 

3.6 
5.6 



B- 



1.1/ j^- 
0.9 V^O 

1.5^0 



j+0.56 
-0.57' 



4.23] 

1.8 ± 0.2(Ff- 
l.l±0.1(Ff- 

1.8 ±0.2 (Ff- 

5.1 ±1.6, 
3.0 ± 0.3(Ff~ 

1.9 ± 0.3(1^- 
3.0 ± 0.3(Ff- 



+o 
-l 

+0 
-0. 
+0. 

-1 



■ K ) 
■ K ) 
■ K ) 

■ K ) 
■ K ) 

>K 



?(t) iloK 
?(7) ±f:l(^ s 
Id) +:iS 



-2.6 (Xa), 
-1.9UAJ, 



-1.9' 



m 



+0.7/ A* 
-0.6 1^0 

± 0A{A*~ 



-0.6 ^ 



±0.5(A 3 " 
±0.2(^ 3 " 
) ±0.5(^" 



-2.6 



Data, 
SM, 

Case — B, 
Case — C, 



(28) 



p ) tli(m s 



\ +2.3 
I -0.4 
3.0 
0.8 



(xa) 
(Xa), 



+1.9, 
-1.1* 
+1.7, 
-0.9 1 



m. 



tit(XA), 

+ -IKxa), 



i^s) tiiixA) 



Data, 
SM, 
Case • 
Case • 

Data, 
SM, 
Case • 
Case • 



,(29) 



B. 
C. 



,(30) 



Here the new dominant error source is the uncertainty of the mass m s , 85 MeV < m s < 125 
MeV. In Fig. IH we show the 7 dependence of the branching ratios for the three measured 
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FIG. 4: The same as FigQbut for B — * Kp decays. 



B — > i^p decays respectively. 

From Eqs. (|28H3T)j) and Fig. HJ one can see that the central values of the theoretical 
predictions in the SM and the mSUGRA model of case C are almost identical and about 
half of the measured values. In case B, unfortunately, the SUSY contribution produce a 
further forty percent variation in the "wrong" direction. But when we include the large 
errors induced by the uncertainties of xa, and form factors, the theoretical predictions 
in both the SM and the mSUGRA model (both Case B and Case C) can be consistent 
with the data. 

For B — ► Kuj decays, the situation is very similar to the B — ► Kp decays. From 
Eqs. ()3HI32|l and Fig. one can see that the central values of theoretical predictions in 
both the SM and the mSUGRA model are much smaller than the measured values. And 
the SUSY corrections in Case B even lead to a further reduction of the branching ratios. 
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FIG. 5: The same as Fig^but for B — ► Klo decays. 



The data and the theoretical predictions (in units of 10 6 ) are 
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Here one can see that (a) the central values of the theoretical predictions in both the SM 
and mSUGRA model are all smaller than the data; and (b) the dominant error still comes 
from the uncertainty of the annihilation contributions and is very large in size. 



4. B -> K m r], B 



V P 



For the decay channels involving rf^ meson, the dynamics is rather complex and has 
been studied by many authors, for example, in Refs. [l^, E3|. Here we didn't consider 

the additional form-factor type contribution to the flavor-singlet coefficients a^Mi^]) ( 
see Ref. [l3| ) since it has large uncertainties. The data and the theoretical predictions 
(in units of 10~ 6 ) for the three measured decay modes in both the SM and the mSUGRA 
model are 
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FIG. 6: The same as FigQbut for B — > K*rj and B —* p t\ decays. 



Br(B u -> K*°7j) 



Br(B — > p rj) 



18.7 ± 1.7 

6.6±0.8(A^ K 
[G.ltlHA^n ±0.1( 7 )^(m s )^( XA ) 



±0.1( T )^(m s ) ^(xa), 
) ±0.1( 7 )^(m s )^( XA ; 



8.6 



+1.9 
1.8' 



6.3 ±0.2(» ±0.1(7) ±0.1(m s ) ±0.3(xa), 
6.4+^0) +^( 7 ) ±0.1(m s ) ±0.3( X a), 
6.3 ±0.20) ±0.1(7) ±0.1(m s ) ±0.3(xa) ; 



DclLcl, 

SM, 
, Case — B, 
Ccise C ; 

SM, 

Case — B, 
Case — C. 



(34) 



(35) 



Clearly the uncertainties of xa arid m s are the dominant error sources here. 
By comparing the numerical results with the measured values, one finds that 

• The SUSY contributions to these three measured decays are always small for both 

*o 

case B and C. For B~ — > K*~r\ and B — > K rj decays, the central values of 
the theoretical predictions in both models are only around 30% of the measured 
values, and the data can be explained only when we take the large theoretical and 
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experimental errors into account. For B~ — > i^~r] and B K r\ decays, however, 
the theoretical predictions in the SM and mSUGRA model as given in Ref. [26] 
agree well with the data. 

• For B~ — > p~rj decay, a b — > d transition process, the SUSY contribution is very 
small, and the theoretical predictions in the SM and mSUGRA model are consistent 
with the data within one standard deviation. 

• As illustrated in Fig. El the branching ratios of all the three channels have a weak 
dependence on the angle 7. For B —>■ K*r] decays, the theoretical errors are large in 
size and induced dominantly by the uncertainty of the parameter \a and the mass 
fn s . For B~ — ► p~i] decay, however, the theoretical errors are relatively small. 



V. BRANCHING RATIOS FOR B — > VV DECAYS 



For B — ► VV decay modes, one generally should evaluate three amplitudes with dif- 
ferent helicity since they all can make contributions and do not interfere with each other. 
In terms of the helicity matrix elements 

H x =<V x {X)V 2 {X)\H eff \B)>, A = 0,±1, (36) 
the branching ratios of B — > VV decays can be written as 

Br(B -> VM) = r B J^L [|tf Q |2 + |# +1 |2 + ^2] (3?) 

The three independent helicity amplitudes Hq, H + \ and H-i can be expressed by three 
invariant amplitudes a\,b\,c\ defined by the decomposition 



(3* 



Here e fJl (r] l/ ), pi^ and Mi^ are the polarization vector, four momentum and masses of V\p, 
respectively, while p = p\ + p 2 is the four-momentum of B meson. The helicity elements 
H\ can be further simplified as 

H±i = a±i ± c±iV x 2 — 1, Hq = —a x — b (x 2 — l) 
M 2 — M 2 — M 2 

x = B 1 2 (39) 
2M l M 2 y 1 



A. Helicity amplitude of B — > p + p , an example 



Now we take the decay B° — > p ± p T as an example to show the ways of decomposition. 
With the QCD factorization approach, the decay amplitude of B° — > p ± p T decay reads 

.Gf 

tin I ~ , 1 1 in , , -\- 1 1) i -! V 'In, , 



A\B -> p + p- 



-i—f p m p 



(e+ ■ e_)(m B + m p )Af^ p (m 2 ) 



2A 2 B ^(m 2 J 



v/2 

-Pfl)(e- "Ps, 



^ V OL (3 



V B ^ p (m 



TUB + rrir 



(40) 
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By comparing this expression of decay amplitude A\ with Eq. (|38j) . one can find the 
coefficients a x , b\ and c\, 



ax(B- 


-P 




b x (B- 




>-) 


cx(B- 


-P 





r* A 



-^J=f p rn p (m B + rr 
V2G F f p m 3 p A2 ' (mp V ^K>i - V tb V td [a x + <]}, 

Tfl q ~\~ /Tip 

V B ^ p (m 2 ) 

' ' {Vu^A- V tb V td [a x + a x }}, 



r * I „ A i A 



where the factorized coefficients a x , a\ and a A can be written as 



a x = d 



C 2 a,C 



+ 



^ ^ [/* + /*] 



iV c 4vr iV c 

C3 a s Cp \ j.\ , r \ -| , oc s Cp 



^G x (s u ) + ^G A (s c 
Vt v t 



+C 3 [G x (s q ) + G x (s b )} + (C 4 + C 6 ) 



G x (s, - 



(41) 



(42) 



-C7 9 [e g G A ( S ,)+ 



e b G x (s b )} + -(C 8 + C 10 )$> q 



(43) 
(44) 



In the above equations, the vertex corrections f x and the hard spectator scattering 
contributions f X j are given by 



fx 
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•12 In — - 18 + / dM< 2 (w) (3- — — In // - 3/tt 
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(45) 
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du- 
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In u — 3in I , (46) 



(47) 



47T 2 2if B f^J V2 m v , 
N c m B h± 

x / du\gf V \u)- 



-(l=Fl) / ^ 



1 /v*^ 



+ 



iV, 



c 



b' l ± 



d^- 



/(a)Vi, 



x d»d» aS-)»±£k_w 9 <«>» ± 



gl [u) \ u + v 

-9 ' 



(48) 



here f = 1 — v and the light-cone distribution amplitudes (LCDAS) $Y(u), §^_(u), g± V , 
g^ V can be found in Ref. (HI ]. 
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The functions describing the contributions of the QCD penguin-type diagrams , th e 
dipole operator Sg and the electro- weak penguin- type diagrams in eqs. (|42MH) are |4l| 



(49) 



G°(s) = + 4 / du<S>f(u)g(u,s), 

3 3 m b J 11 

99 rl i(a)V'2 / \ 

G±(s) = Z --Un^- + 2 d^fffWi^Mv), (50) 
3 3 m b J - 4 



G° 9 
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i 2^ 2 (u) 
du 



o 



1 -u 



-J\u ^ )V \u) 
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l-u 



- na)V2 ( \ 

- (v)v 2 1 s , u 9± W 
-ug\_' (u) + 



{a)V2 ( \ 

9± W 



^G\s u ) + ^G\s c 
v t v t 



C 2 + ^ 



N, 



c 



with the function g(u, s) defined as 



g{u,s) 



dx xx In (s — xxu — ie). 



(51) 
(52) 



(53) 
(54) 



(55) 



From Eq. (j39J) . we can finally get the helicity elements H±i and H . One should note 
that the coefficients a, b and c are independent of the helicity A in the naive factorization 
approach. In QCD factorization approach, however, the coefficients ax, b\ and c\ depend 
on the choice of A (A = 0, ±1). 



B. Numerical results 

Using the above formulae it is straightforward to calculate the branching ratios of 
B — > VV decays. In Table I1VI and |VJ we show the theoretical predictions for the CP- 
averaged branching ratios of the nineteen B — > VV decays in both the SM and the 
mSUGRA model, assuming = mf ) /2,mf, and 2m\ ) , respectively. 

In numerical calculations, we do not consider the annihilation contributions and only 
give Br-f . The reasons for this can be found in the last paragraph of Sec. IIHI Of course, 
the annihilation contribution may play an important role to explain the large transverse 
polarization of the B — > K*°p + and (f)K* + decays, as pointed out in a a recent paper 
[44J , where the author has tried to include the contribution induced by a QCD penguin 
annihilation graph to explain the observed longitudinal polarization, fliB —* <pK*) « 
50%, in the framework of the SM. More studies are needed to understand this problem 
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TABLE IV: Numerical predictions in the SM and mSUGRA model (both case B and case C) 
for CP-averaged branching ratios in units of 10~ 6 for b — > d transition processes of B — > VV 
decays, here the annihilation contributions are not included. 
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clearly. But it is beyond the scope of this paper, and we here just calculate the possible 
new physics contributions to the CP-averaged branching ratios of B — > VV decays in the 
mSUGRA model, and show the general pattern of such contributions. 

From the numerical results as given in Table IIVI and one can see that the new 
physics contributions to the branching ratios of B — > VV decays are generally small for 
the SM-like case-C, but can be rather large for the case-B, especially to those penguin 
dominated B — > K*(p, u,<p) ( b — > s transition ) processes. 

Among the nineteen B — ► VV decays, only seven of them have been well measured, 
while upper limits are available for the remaining twelve decay modes. For the world aver- 
age of the experimental measurements, one can see HFAG home page 4^| and references 
therein. In the following subsections we focus on the seven measured decay channels. 
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TABLE V: Same as Table HVl but for b — ► s transition processes of B — > VV decays. 
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1. B»^p + p-,B 



These three decays are tree-dominated decay channels. The data and the theoretical 
predictions of the branching ratios (in units of 10~ 6 ) in the SM and mSUGRA model are 
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(57) 
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Here the dominant errors also come from the uncertainties of the form factors and the 
angle 7 

Fig. [7| shows the 7 dependence of the branching ratios for the three decays in the SM 
and the mSUGRA model. From this figure and the numerical results as given in Eqs.([561 
I58|) . we can easily find that the SUSY contributions to these tree-dominated decays are 
very small. The largest error here still comes from the uncertainty of the relevant form 
factors, and the theoretical predictions in both the SM and mSUGRA model are consistent 
with the data within one standard deviation. 
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FIG. 7: The same as Figdbut for B° -> p + p~ and 5" 



~(p ,cu) decays. 



2. 



B- 



K*-p°, K*°p- and B 



K* 



Since these four decays are penguin dominated decay modes, the SUSY contributions 
may be significant. The world averages, the theoretical predictions (all in units of 10~ 6 ) 
in the SM and the mSUGRA model are 
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FIG. 8: The same as Figfflbut for B~ -> K*-p°, K* p~ and B -» K*<f> decays. 
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( 4.6±i:|0*) ±H(A B ^ K *) ±0.02(7), Case-C. 

Here the major errors are induced by the uncertainties of the input parameters /i, A B ^ P , 
A B ^ K and the angle 7. 

From the numerical results as given in Eqs. (j59H62j) and Fig. |S1 one can see that 

• For these four channels, the central values of the SM predictions are less than half 
of the measured values. The SUSY contributions are negligibly small for the case 
C, but can provide a 70% enhancements for case C. 
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• The dominant errors still come from the uncertainty of the form factors and the 
renormalization scale \i. Except for the decay B~ — > K*~ p°, other three decays are 
almost independent of the angle 7. The reason is that the latter three decays are 
pure penguin processes and have no terms proportional to V t d or V U b. 

VI. SUMMARY 

In this paper, we calculated the new physics contributions to the CP-averaged branch- 
ing ratios of the thirty nine B — > PV and nineteen B — > VV decay channels in the 
mSUGRA model by employing the QCD factorization approach. 

In Sec. |nj a brief review about the mSUGRA model and the SUSY corrections to the 
Wilson coefficients was given. In Sec. IIH| we presented a short discussion about the QCD 
factorization approach which is commonly used to evaluate the hadronic matrix elements 
< MiM 2 \Oi\B >. In Sec. II VI and Sec. |V| we calculated the branching ratios of B — > PV 
and B — > VV decays in the SM and the mSUGRA model, and made phenomenological 
analysis for those well-measured decay modes. 

In the mSUGRA model, the Wilson coefficient CV 7 (mh) can be either SM-like or sign- 
flipped comparing with that in the SM. The data of B — > A s 7 can provide a strong 
constraint on the size of CV 7 (mb), but not its sign. The latest measurements of B — > 
X s l + l~ decay prefer a negative ( SM-like) CV 7 (mb). In this paper we choose three typical 
sets of the mSUGRA input parameters (m , mi , A , tan/3, Sign(p)) in which CV 7 (m{,) can 
be either SM-like ( the case A and C) or have a flipped-sign (the case B). 

As expected, the SUSY contributions to all channels considered here are very small 
for both the SM-like case A and C. For case B, however, the SUSY contributions can be 
significant for those penguin-dominated decays, and for this case we found that 

• For those tree-dominated decays or the decay channels having the penguin contri- 
bution only coming from a^, a™ or having only weak annihilation contribution, the 
SUSY contributions in case B are also very small. 

• For those QCD penguin-dominated decay modes, the SUSY contributions can in- 
terfere with the corresponding SM parts constructively or destructively, and conse- 
quently can provide an enhancement about 30% ~ 260% to the branching ratios of 
the decays B — > nK*,K(j),K*(j),K*p etc, or a reduction about 30% ~ 80% to the 
branching ratios of the decays B — > Kp, Kuj etc. 

• For B — > K*{ji,p) and B — > (K,K*)(f) decays, the central values of the SM predic- 
tions for branching ratios are generally less than half of the measured values. But 
the mSUGRA predictions can become consistent with the data within one standard 
deviation due to the inclusion of the significant SUSY enhancements. 

• For B — > Kp and B — > Kuj decays, the SUSY contributions interfere with their SM 
parts destructively, more investigations are needed to cover the gap between the 
theoretical predictions and the data. 

• For most decay channels, the error induced by the uncertainty of the annihilation 
contributions is dominant and large in size. The SM predictions, consequently, have 
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large theoretical errors when the QCD factorization approach is employed. The 
large SUSY contributions in the case B , unfortunately, may be masked by the large 
theoretical errors dominated by the uncertainty from our ignorance of calculating 
the annihilation contributions in the QCD factorization approach. Hence we expect 
the theory will go further to make the annihilation contributions calculated reliably 
and only at that time can one distinguish the new physics contributions from the 
theoretical errors. 

In short, only in case B, can the SUSY contributions be significant to the branching 
ratios of the studied B meson decays. If the sign of CV 7 (m&) can be strictly determined 
to be SM-like finally, as claimed firstly in Ref. [37], the parameter space of the mSUGRA 
model will be strongly restricted, and the branching ratios of the two-body charmless B 
decays can hardly be affected. 
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APPENDIX A: INPUT PARAMETERS 

In this appendix we present the relevant input parameters being used in our numerical 
calculations. 

• Decay constants. The decay constants are defined by following current matrix 
elements 1 451: 



(P(Q)\qwq\0) = (V{q,e)\q^q\0) = fvm v e*». (Al) 

And the scale-dependent transverse decay constant in Eq.()13|) is defined as[l3| 

(V(q,e*)\qa^q\0) = /v(<7„< - <fcej) , (A2) 

where e* denotes the polarization vector of the vector meson V. These decay con- 
stants are nonperturbative parameters. And they can be extracted from the exper- 
imental data or estimated with well-founded theories, such as QCD sum rules, etc. 
The decays constants we used here, in units of MeV, are 



/tt I k fs 


fp fK* fu) fcf) 


f-L f-L f-L f-L 

J p J K* J uj J q> 


131 160 200 


209 218 187 221 


150 175 150 175 



As to the decay constants related to rj and r\ , we shall take the convention in Ref. 

^\qipimW\p)) = -ifl^Pn ( A3 ) 
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the quantities /" (/) and in the two-angle mixing formalism are 

= A S i n e 8 + A cos ^ , /- = _ 2 A sin # 8 + _| cos ^ (A4 ) 

f u = JjL cos e 8 _ljL gin ^ ^ f s = -2 A cos ^ 8 _ J| sin # ( A5 ) 

with f 8 = 1.28f n , f = 1.2f v , 6 = -9.1° and 6 8 = -22.2°. In order to incorpo- 
rate the charm-loop diagram contribution into rf^ involved decays, two new decay 
constants should be given, namely /° and f c , . Following Ref. |l2j, we take values 
L c ~ -lMeV and f c , ~ -3MeV. 

Form factors. The form factors are in nature nonperturbative quantities, extracted 
usually from the experimental measurements. In Ref. 45||, the form factors were 
defined by current matrix elements 



(P(q)\?f(l-y B )q\B) 



P B + <? - -^-/-V 



m B — rrii 



k 2 



(A6) 



(F(g,e)|^(l-7 5 )g|B> = e^e w p a B q 



m B + my 



A' 2 



+ ie*(m b + my)Ai(/c 2 ) - z 



mft + my 



(pB + g)„A 2 (fc 2 



(A7) 



where k = p B — q. We neglect corrections to the decay amplitudes quadratic in the 
light meson masses, so that all form factors are evaluated at k 2 = 0. At the poles 
k 2 = 0, we have 



F (0) = Fx(0), A (0)=A 3 (0), 
2m v A 3 (0) = (m B +m v )A 1 (0) - (m B -m v )A 2 (0). 



(A8) 
(A9) 



In our calculations, we use the form factors as given in Refs. [13j,|46[. They are 



B 



7T 



B -> K 



B^p 



B ^ K* B -> 



V 



0.28 ±0.05 0.34 ±0.05 



0.37 ±0.06 0.45 ±0.07 0.33 ± 0.05 
0.242 ± 0.02 0.292 ± 0.06 0.219 ± 0.02 
0.221 ± 0.02 0.259 ± 0.06 0.198 ± 0.02 
0.323 ± 0.03 0.411 ± 0.08 0.293 ± 0.03 



For the form factors of B — > rj and B — > r\ transitions, we use 



r 0,l ~~ r 0,l 



. cos sin 9n 



V6 \/3 



F W =F B f . (! m|s + ^. ) (A10) 
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The parameters Xh,a- When calculating the contribution from hard spectator scat- 
tering and annihilation diagrams, the end point divergence will appear. To treat 



11 



And 



the divergence, one can parameterize them by \h and \a respectively [11 

In numerical calculations, we set pn = Pa — as the default input values 
we also scan over pA € [0, 1] and 4>a € [—180°, 180°] to estimate the theoretical error 
induced by the uncertainty of the annihilation contributions. 



In this paper, we use the same CKM angles, quark masses, meson masses and the 
B meson lifetimes as used in Ref. 12611 . 
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